Small World Patterns in Food Webs 
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The analysis of some species-rich, well-defined food webs shows that they display the so called 
Small World behavior shared by a number of disparate complex systems. The three systems analysed 
(Ythan estuary web, Silwood web and the Little Rock lake web) have different levels of taxonomic 
resolution, but all of them involve high clustering and short path lengths between species. Addition- 
ally, the distribution of connections P(k) is skewed in all the webs analysed and shows a power-law 
behavior P(k) oc fc~ 7 in two cases (with 7 « 1). These features suggest that communities might be 
self-organized in such a way that high homeostasis to perturbations (with short transient times to 
recovery) would be at work. The consequences for ecological theory are outlined. 



I. INTRODUCTION 



The understanding of food web patterns in complex 
ecologies is a central issue in ecology. As Stuart Pimm 
claims, food webs are the road-maps through Darwin's 
entangled bank (Pimm et al, 1991). Beyond the spe- 
cific features of currently described food webs (such as 
the concrete species composition) a number of regulari- 
ties have been observed suggestive of fundamental laws 
of ecosystem organization (see for example Warren, 1994; 
Polis and Strong, 1996 and references therein). 

One of the most successful approaches to real food 
webs deals with the graph patterns associated with the 
set of links between species. Specifically, a food web can 
be described in terms of a graph G(V, E) consisting of a 
finite set V of vertices (nodes, species) and a finite set E 
of edges such that each edge e is associated with a pair of 
vertices v and w. In spite that ecological graphs are di- 
rected (i. e. links go from one species to another and the 
reverse connection will be typically different in strength 
and sign) most theoretical studies deal with non-directed 
graphs. In other words, two species appear connected 
if they share a common edge, irrespective of its particu- 
lar properties. This approach has been rather successful, 
since many reported regularities form field data are fully 
recovered from non-directed graph arguments of different 
nature (Cohen et al., 1990; Pimm et al., 1991; Williams 
and Martinez, 2000; see also Albert et al. 2000). This ob- 
servation is probably a consequence of fundamental laws 
constraining the basic topological arrangements of eco- 
logical networks that go beyond the specific rules of dy- 
namical interaction (Brown, 1994; Sole and Bascomptc, 
2001). 

Recently, new theoretical approaches to graph com- 
plexity in nature have emerged. Two main results of 
these studies are: (a) the widespread presence of the 
so called small world (SW) pattern of some economic, 
technologic and biological networks (Watts and Strogatz, 
1998; Watts, 1999; Adamic, 1999) and (b) the pres- 



ence of scale-free (SF) distributions P(k) of connections 
(Barabasi and Albert, 1999; Jeong et al., 2000; Wagner 
and Fell, 2000). Specifically, in some cases the number 
of nodes with k links follows a power law distribution 
P(k) ~ k~ 7 where most units are weakly connected and 
a few are highly linked to other nodes. This is however 
a non-universal feature: the graph of the neural network 
of C. elegans, for example, displays an exponential dis- 
tribution of connections (Amaral et al., 2000). 

The SW pattern can be detected from the analysis of 
two basic statistical properties: the so called clustering 
coefficient C v and the path length L. Let us consider the 
set of links = 1, S) where S is the number of 

species and let's assume that £y = 1 if a link exist and 
zero otherwise. Species will be labeled as Sj(i = 1, S). 
If a food web is considered, an additional parameter is 
the community matrix connectance C. Let us consider a 
given species (the i-th one) and the set of nearest neigh- 
bors Ti = {s.^ij = 1}. Here we can calculate the cluster- 
ing coefficient for this species as the number of connec- 
tions between the species belonging to F^. By defining 
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and thus c v (i) — £i/(CS(CS — l)/2) so that the cluster- 
ing coefficient is the average over all species: 



C v — — ^ c v (i) 
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and measures the average fraction of pairs of neighbors 
of a node that are also neighbors of each other. 

The second measure is easily defined. Given two ar- 
bitrary species Sj and Sj, let L m i„(i,j) be the minimum 
path length connecting these two species. The average 
path length L will be: 
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Small World graphs are highly clustered but the min- 
imum path length between any two randomly chosen 
nodes in the graph is short. By comparison, random 
graphs (where nodes are randomly connected with some 
probability) are not clustered and have short L (Watts, 
1999). At the other extreme, regular lattices are typically 
clustered and have long distances. It has been shown, 
however, that a regular lattice can be transformed in to 
a SW if a small fraction of nodes are rewired to randomly 
chosen nodes. Thus a small degree of disorder makes the 
lattice to have short paths (as in the random case) but 
is still mostly regular (Watts and Strogatz, 1998). 

In a graph with an average of < k > links per node, it 
can be shown for random graphs that C™ nd ~< k > /N 
and L rand will be short. For large networks, a SW is 
present if L rand is larger, but close to L and C v > C r v and . 
When networks composed by a small number of units 
are analysed, the second condition is often replaced by 
C v > C r v and . This is the case for the SW pattern in the 
neural network of the nematode worm C. elegans (ce) 
which has 282 neurons and C£ e = 0.28 and L ce = 2.65, 
to be compared with the corresponding random graphs: 
C vand = o_Q 5 and L rand = 2 2 5. For the metabolic net- 
work of E. coli, Wagner and Fell obtained similar results, 
but with a much larger difference in the clustering. 

The consequences of the SW and SF paterns are far 
from trivial and can be of great importance in recogniz- 
ing evolutionary paths, the origins of homeostatic stabil- 
ity and the sensibility to perturbations in biological net- 
works. Watts and Strogatz discuss some of these ideas 
in their seminal work, suggesting that a SW architecture 
would play a relevant role in enhancing synchronization 
in the visual cortex. For metabolic networks, Wagner and 
Fell suggest that it might play a relevant role in allowing 
metabolism to react rapidly to perturbations thus dis- 
playing a very high homeostasis. What about ecological 
networks? Do food webs display this type of topological 
properties? 

In this paper we explore this question by analysing the 
statistical properties of a set of large (S > 100) ecological 
networks where a fine taxonomic resolution is available. 
For these webs, the SW property is shown to be present 
and a skewed distribution of links is also shown to be 
involved in most cases. The presence of these properties 
will be shown to have deep implications for theoretical 
and applied issues at different levels. 



II. ECOLOGICAL SMALL WORLDS 



tures with some standard graph patterns such as Cayley 
trees (Watts, 1999). For the constant connectance model 
(Martinez, 1991a) no such ranking is introduced and a 
trophic species preys on any other trophic species (in- 
cluding itself). However, ecological nets display both 
types of properties and some studies suggest that they 
also display clustering. 

Evidence for clustering in food webs has been com- 
piled in recent years (see Solow and Beet, 1998). A num- 
ber of authors have shown that food webs are typically 
non-random but the effects such topological properties 
on stability are far from clear (see Solow et al., 1999, for 
a recent discussion). 
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FIG. 1. (a) clustering coefficient for the four analysed food 
webs (see text). Here the dark bars are real webs and white 
bars correspond to randomly generated webs with the same 
average number of links per species (averaged over 200 gen- 
erated samples). We can see that in all cases the clustering 
is clearly larger than random; (b) Characteristic path length 
L with random values. Except for Little Rock, the difference 
between the random and real case is very small. For the two 
Ythan webs, the improvement in resolution slightly modified 
the observed relations towards a better defined SW pattern. 



It is interesting to see that in fact most static graph 
models of food webs involve features ranging from purely 
random nets (constant connectivity networks, see Mar- 
tinez, 1991a) to purely hierarchical models (the cascade 
model, see Cohen et al., 1990). Under the second class, a 
trophic species can only prey on a trophic species of lesser 
rank leading to a hierarchical structure sharing some fea- 



The possible importance for clustering on the stability 
of ecological networks was already recognized from early 
statistical approaches (May, 1974; Pimm, 1982). Given 
the potential implications of the SW phenomenon for net- 
work homeostasis, the presence of such properties in real 
food webs would be of great interest in this discussion. 
We have found that, in fact, rich-species food webs with a 
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good taxonomic resolution display the properties of small 
world behavior. In figure (I) we show the obtained results 
for the four large food webs analyzed, which are known 
in the literature because they are the most detailed large 
community food webs available. These are: 




Number of links k Number of links k 



FIG. 2. Frequency distributions of links P(k) for the anal- 
ysed webs. Left column: actual data, which fit quite well a 
power law for the Ythan and Silwood cases (the exponen- 
tial fit is also shown as a dotted line). The improvement in 
taxonomic resolution for Ythan (2) only modifies the expo- 
nent slightly. Right column: Expected distributions from the 
random graph approximation, predicting Poissonian behav- 
ior. None of the analysed food webs can be approached by 
the random case. 

1. Ythan estuary: this web is from a freshwater- 
marine interface. The number of species reported 
from two different studies are S = 93 (Hall and 
Raffaelli, 1991) and S = 134 (Huxham ct al., 1996). 
These studies differ in the taxonomic detail incor- 
porated and we use them in order to check the sen- 
sitivity in relation with the measurement of SW and 
scaling properties. This is the largest documented 
web in the UK and the average number of links per 
species are < k >= 8.75 and < k >= 8.70, respec- 
tively. The second web was expanded by adding 42 
metazoan parasite species but their effect on < k > 
was weak. For this web most (88%) nodes corre- 
spond to real species. Some nodes involve lower 
resolution at different levels (Nematodes, Acarina 
or Brown algae, for example, are lumped together). 

2. Silwood park web (Memmott et al., 2000): this is 



a sub-web obtained from a field site of 97 hectares 
in size. This is a very detailed community cen- 
tred on the Scotch broom Cytisus scoparius, in- 
volving S — 154 species: one plant, 19 herbi- 
vores, 5 omnivores, 66 parasitoids, 60 predators 
and three pathogens. The average connectivity is 
< k >= 4.75. All nodes of this (sub-) network are 
species (except node 122: immature spiders). 

3. Little Rock lake (Martinez, 1991): a web from a 
freshwater habitat, it includes S = 182 species. 
The Little Rock lake is a small lake with an area of 
18 hectares. Here < k >= 26.05. For this web, only 
31% of nodes are associated to species. Most are 
genera-level nodes (63%) and the rest correspond 
to higher taxa (Bivalvia, Hirudinea, etc). 

We compared the observed values of C v and L with the 
expected ones from the randomized webs with the same 
number of (total) links. As expected from a SW pattern, 
we see that all have very similar (and very short) dis- 
tances L and a clustering coeficient from C Ythl / 'C r v and = 
2.33 to C^ ilw /C r v and = 5.0 (numerically close to the C. 
elegans data) i. e. clearly larger than expected from ran- 
dom webs. For the average path length, we found that 
the values are very close except for the Little Rock case, 
were notable differences are present. 

An additional and surprising feature of the observed 
distributions of connections P(k) is that they show 
strongly non-Poissonian behavior (wich would be ex- 
pected from random wiring). The classic result by Erdos 
and Renyi (ER) on random graphs (see Bollobas, 1985) 
with S nodes shows that the probability that a ver- 
tex has k edges follows a Poisson distribution P(k) = 
e~ <k> X k jk\. Here 

<k>=(^ S ~ rS jP k (l-P r ) s - k - 1 (2) 

where P r is the probability that two nodes are connected. 
The observed distributions for Ythan webs and the Sil- 
wood example are shown in figure 2 (left column). In 
both webs, the distribution was found to be strongly 
skewed with a good agreement with a power law fit. The 
observed values for a regression of the log-transformed 
data are: 

j Ythl = 0.94 ± 0.06 (r 2 = 0.79,p < 0.01) 

j Yth2 = 1.04 ±0.05 (r 2 = 0.83,p < 0.01) 

7 Silw = 1.13 ±0.06 (r 2 = 0.79, p < 0.01) 

For comparison, the expected connectivity distribu- 
tions from a random, Poissonian net are also shown in 
Figure 2 (right column). Here the average number of 
links per species in each web was used. 
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The Little Rock web (figure 3) shows a fluctuating dis- 
tribution of connections with a very high variance, but 
no obvious recognizable standard shape, although it is 
also skewed, with a prominent peak of 21 species (involv- 
ing a coarse-grained cluster of cyanobacteria and green 
algae predated by pelagic cladocera, copepods and ro- 
tifers) with k = 26 links. These anomalies are likely to 
due to the taxonomic resolution used (mostly at genera 
or higher levels). 

The exponent derived from our study is rather dif- 
ferent from the estimations obtained for Internet or the 
metabolic graphs suggesting that, in spite that the SW 
pattern together with the scaling law might be common 
to many different types of networks, the ecological webs 
belong to some different class of complex phenomena 
(Sole et al., 1996) to be defined in future theoretical stud- 
ies. 
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FIG. 3. Little Rock lake link distribution. Here P(k) is 
clearly skewed but strong fluctuations are observable as scat- 
tered peaks, in particular a large peak at k = 26 involving 
a large cluster of interacting species (see text). Other peaks 
are indicated, revealing strong deviations from a monotonous 
power-law tail. 



III. DISCUSSION 

In this paper we have presented evidence for SW pat- 
terns in ecological networks. The aim of our study was to 
see to what extent such evidence was strong enough by 
analysing a small number of large, taxonomically well- 
defined food webs from the ecological literature. The 
webs analysed have different origins but all of them share 
the SW pattern in their graphs and are clearly away from 
their random counterparts. As far as we know, this is the 
first reported evidence for such property in ecology. It re- 
veals, by means of a new measure, the existence of clus- 
tering in food webs associated to short path lengths and 
opens a number of questions not previously formulated 



in previous theoretical studies. 

Although some food web analyses found evidence for 
the existence of compartments (Solow and Beet 1998), 
such species clusters were established in terms of the 
trophic similarity between pairs of species (see also Pimm 
and Lawton 1980, for a similar measure). In contrast, 
the aim of the clustering coefficient as a measure of com- 
partmentalization is to find out to what extent groups of 
species are more connected internally than they are with 
another species or groups of species. In other words, this 
new measure go beyond by looking not only for resource- 
use-guilds (see Wilson, 1999) but also for sub-webs with 
higher internal connectance than the rest of the commu- 
nity. 

The use of the clustering coefficient is supported by the 
fact that groups of species with diets of similar species 
do not occur, especially in species-rich systems, such as 
those analysed in this article (Winemillcr 1990; Polis 
and Strong 1996). In this respect, this measure of com- 
partmentalization minimize the importance of discrete 
and homogeneous trophic-level dynamics in complex food 
webs, an issue widely considered as an artifact by many 
ecologists (see Polis and Strong, 1996, for a review). The 
clustering coefficient does not consider the trophic posi- 
tion of the species in the food web, focusing only in the 
presence or absence of the interactions. 

Two of the webs analysed correspond to the same sys- 
tem: the Ythan estuary food web. These two webs in- 
volve different taxonomic resolution but the exponents 
and indices remained basically the same for the two cases. 
This agreement suggests that the SW property and the 
scaling law are rather robust and that no completely de- 
tailed description of the species-level graph needs to be 
defined in order to detect these properties. For the best 
defined web (Silwood), where all nodes belong to species, 
we obtain the best evidence of SW, suggesting that less- 
detailed webs introduce a coarse-grained resolution that 
might hidden a higher degree of clustering (this seems to 
be confirmed by the Ythan example). The worst result is 
given by Little Rock lake, where the resolution is mainly 
at the genus level and strong deviations from the power 
law are observable. 

The SW phenomenon appears to be associated with 
highly non-random, fat-tailed distributions of links. The 
plots of P(k) show that all food webs (and others not 
explored here due to their lower taxonomic resolution) 
involve strong deviations from the Poissonian behavior 
and confirm the presence of clustering in ecological inter- 
actions. These properties reveal that some new features 
of ecological graphs have to be considered by future mod- 
els in order to explain their origin and meaning. 

In our system, the presence of such a web pattern 
might have important consequences for stability and 
resilience. Species removal (Brown and Heske, 1990; 
Brown, 1994) or the introduction of sharp shocks in 
population abundances due to transient perturbations 
might be strongly damped through the rapid transfer 
of information characteristic of these webs (Watts and 
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Strogatz, 1998). In this sense, classic studies on ecosys- 
tem resilience and resistance might benefit from consid- 
ering these topological constraints. Most studies sug- 
gest that the characteristic time of recovery from short, 
sharp shocks (as it would be the case for a pulse of nu- 
trients) depends on the lenght of the food chain. Some 
enlightening theoretical studies based on the analysis of 
small-species models show this trend (Pimm and Law- 
ton, 1977) where some simple topological arrangements 
between species are predefined. But perhaps some other 
regularities will emerge from our results when rich-species 
models are taken into account. 

A different question emerging from our study is how 
such food webs are constructed through community as- 
sembly rules (Drake, 1990a, 1990b). Since membership in 
a community is usually influenced by interactions among 
species, the dynamical process leading to a SW pattern 
of connections might eventually allow to understand how 
our reported results are explained. A possible dynam- 
ical process underlying the presence of scaling in food 
webs is preferential attachemcnt (Barabasi and Albert, 
1999). This process of assembly consists in the addittion 
of new nodes which arc preferentially linked to already 
highly-connected nodes. This process has been argued 
as the explanatory mechanism of the organization of dif- 
ferent complex networks, including technological ones or 
the large set of metabolic networks analysed by Jeong 
and collaborators (2000). But for ecological networks 
(represented as food webs) is far from trivial to decide 
if preferential attachement is an appropriate, general as- 
sembly rule. The ubiquity of assembly rules has been 
recently discussed (Weiher and Keddy 1999) almost each 
one corresponding to each of the different ecological com- 
munities studied. 

Even if a community is assembled by random addi- 
tion of species, interactions can lead to non-random webs. 
If network homeostasis is a global constraint to species 
richness, then a SW pattern might be an 'attractor' in 
community assembly (Weiher and Keddy, 1999). Fu- 
ture models will require to incorporate weighted, directed 
links and the scaling behavior of these graphs can have 
some important effects in extinction dynamics (Sole and 
Montoya, in preparation) and might lead to some unex- 
pected results. The SW topology might be related to 
underlying community level constraints and might help 
understanding how biodiversity emerges and persists. 
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